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Abstract 

We present the results of the computation of a twisting type N solution to vacuum 
Einstein equations following an iterative approach. Our results show that the higher 
order terms fail to provide a full exact solution with non- vanishing twist. Nevertheless, 
our fourth-order solution still represents a regular and twisting type N solution. 
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1. Introduction 

There is an interesting unsolved problem of finding solutions of vacuum Einstein 
equations that represent gravitational radiation outside a bounded source. The non- 
twisting solutions of type N, either approximate or exact, have always singularities 
in space at infinity, hence they are not good candidates for this purpose. Also, there 
have been various different efforts to find twisting type N solutions [1, 2] and, so 
far, there is only one exact solution with non- vanishing twist reported by Hauser [1]. 
Mauser's solution also has singularities which extend to spatial infinity. Thus, it fails 
to provide an exact description of the gravitational radiation. 

Recently, Stephani [3] studied a linearised form of twisting type N solutions and 
claimed that either there will always exist singular lines in space for twisting solutions 
or that twisting type N fields do not describe a radiation field outside a bounded 
source. Later on, in [6] an iterative approach has been proposed in an attempt to 
solve this puzzle. 

Using this iterative approach, the linear solution in [3] has been reproduced and 
higher order terms of this iteration process have been described in [6]. In particu- 
lar, the authors have used a first-order solution, flat and non- twisting, to derive an 
approximate non-flat twisting solution. This solution is given by regular functions 
for sufficiently large spheres (i.e. spheres defined by the (C, C) coordinates for large 
values of the affine parameter r). This is an interesting example of twisting solutions, 
where behaviour is contrary to the expectations generated by Stcpliani's paper. 

After all this, one is wondering whether this iterative approach could shed light on 
the general behaviour of solutions describing gravitational radiation from a bounded 
source. In this note, we present the solution obtained by iterating the third-order 
solution in [6] up to higher orders. Unfortunately, nothing can be said about the 
existence of a regular exact solution with a non- vanishing twist parameter. 

2. The vacuum field equations 

Following the notation in [6], we write the metric g in terms of a complex null 
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tetrad, as follows: 



g = Qij^ue'^^e" ^2e^®e'^ + 2e^®e'^ , = e\ and , e"^ real, (1) 

s s 



the overbar indicates complex conjugation. For any vacuum, type-N space-time with 
non- vanishing complex expansion, i.e., where Z = — r42i ^ 0, the null tetrad can 
always be written as [4] 

= -^dC , =idC, -e^ ^du + LdC + LdC , 
PZ ^ ' PZ s s , 

= dr + WdC + WdC-He^ , 
where the metric functions are given by 

= (r - iE) , 2iE = P^(DL - DL) 

W = -Z~^L^u + , D = d^-Ldu (3) 

= -ra«(lnP) + \K, = 2P2 Re[P)(D In P - X„)] 

where the subscripts indicate partial differentiation. Within this tetrad, setting P = 
V^ui the remaining Einstein vacuum field equations take the form 

D[p-^duDDv] = Q , 
DDDDV = DDDDV . 

The curvature and twist are given below and we must require not to be zero 



icW =-^4 = -ZP^dudu{p-^DDV}^0 , 

_ (5) 
2i E =p2 (P)L - DL) ^ . 

The particular form of this (non-zero) curvature component is however not an in- 
variant of this problem, since it may be changed, remaining non-zero, by a choice of 
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transformation to new (allowed) coordinates. It is known that for type N spacetimes 
there are no non-zero invariants of the curvature tensor, which requires us to use the 
gauge invariants of the problem. Therefore, we should use the two (gauge) invariants 
of this problem, the ratios E/r and K/r'^ [5] with the condition that E and K be 
regular functions of their arguments, (see also [7]). 

Fixing the gauge such that 

p = 1 + icC; v = Pu (6) 

and defining 

$ ^L,,+-^^L-LduL (7) 

then the remaining Einstein vacuum field equations (4) can be written as (see [6] for 
details) 

^9„$-l5> = . (8) 

Im I + ^^c^ - '^du^ - Ld^du^ ^ = . (9) 

The N-th order approximation will then be one that approximates the complete so- 
lution, up to that order of iteration, as a sum of all the preceding orders 

AT N 

^^Yl ^^^^ ' ^ ~ ^^^^ (10) 

where L^^^ — and ^^^^ — give the flat space-time. The case = 1 accounts for 
the linearised version in [3] where and L^^) can be obtained from the equations 

d^du^^^^ = and p-2|p2/^(i)| =$(1). (11) 
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At the n-th level of iteration, n > 2 we have 

n-l 

p-2 |p2 ^(n) I ^ Q^^(n) ^ (^^^^^ = $W + J] L^""-^^ d^L^^^ 



^ ^ (12) 

n-l ^ 



whose solutions are obtained from straightforward calculations up to some arbitrary 
functions that can be fixed with the eqn. (9) constraint. Using (5), the curvature 
may be explicitly re-written in terms of the current variables: 

r*4 = p2^1 + iS/r)"^a„2$ . (13) 

the n-th level of the gauge invariants S = ^ T,^'"'^ and K = K^"^^^ at the lowest 
order terms are E^*^-* = and K^^^ = 1. The higher order terms are given by 

e(i) = PHui (14) 
= _2p2Re [d(^duL^^\ (15) 

and for n > 2 

n-l 

= p2|j^ [a^L(")] -Y^l^ \&~^^duL^^^\ } (16) 

n-l 

K^n) ^ _2p2|Re [d^duL^''\ - 5^ Re [L^"-^')^^!^^'^] } (17) 

Next, given a solution at first order one could in principle compute a series that may 
converge or not. 
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3. Fourth-order solution 



Starting with the ansatz, 

lW = 1 and $« = ^(^^) (18) 

and following the iterative approach described above, the authors of [6] computed a 
twisting solution to the third level of iteration. Now, using (12) and solving for L^^) 
and $(4) we have the fourth-order terms given as 

^(2f(^)(.)^) - ^±iflLaa)^(a(i) -^) (19) 

(l+iCO^^^ du ^ C(1 + KC)^ du^ ^ ^''^ 

+ C(i + icc)^ - ^T^ ^ '"^^^^ - T^'"^^^^^ ] 



that satisfies the constraint (9). Now, one can compute E^^^ and K^'^^ from (16) and 
(17), thus 



= iL4^Im (/(4)(«)) _ 2(Re a(i))(Im ^) - 2(Re a(2))(Im ^) 
(l + ICC) ^ " ^ du ' ^ du ' 

" (ITIc?)^'" ^^^^ + ITIc7 ^ d^^^'"^ " ^ 

,Im (7^^) - a(l)^Im (4?)[21n(CC) + i^^(ln(CC)) 



2l 



(1 + iCC)' dx. ^ ^ du^ " 'i + KC 

(21) 



(1 + iCC) ^ ^ 

+ 2a(i)^Re (^)[21n(CC) + ii^(ln(CC))2] 

(22) 

Next, adding all the preceding orders of interation we have the fourth-order ap- 
proximation as follows 



$ ~ - "^ + -')- "^'^"^'^ ^'^^'^ Infm f24) 



%^Im (7M^) - aW^Im (^)[21n(CC) + i|^(ln(CC))^] 



(25) 
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— t^-e(f).4(Re.)(Re0) 

+ (1 + iCC)^ ^ du^^ (1 + KO^ ^ ^ 
+ 4^i^Re(^^^Z!!) 

+ 2a(i)^Re (^)[21n(CC) + i}^(ln(CC))'] 



(26) 



with 



i=i j=2 

In order to avoid the singularities deriving from the logarithmic terms, we can choose 
the first term of the function /2 as /^^^ = iku with k a real constant and leaving /(^^ 
and /*^^^ as arbitrary functions of u. Then, we have a regular fourth-order solution 
given by 



E « M ^l^fg" - 2(Re a(i))] + (/sW) - 2(Re aO(Im ^) (29) 



and 



^,^~^d±Mf^^ (31) 
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with 

At this point we have a non-flat twisting and regular solution of type N. We have 
assumed that there are no non-zero invariants of the curvature tensor and used only 

the gauge invariants. Further calculations indicate that from the ansatz (18) we have 
started with it is not possible to generate well behaved higher order corrections of 
this solution. A new ansatz is needed in order to continue the study of the twisting 
type N solutions using this iterative method. 

We can only conclude that there is not yet exact and regular twisting type N 
solutions. We would need to improved our methods to seek for regular twisting 

type N solutions or to find a consistent proof of the non-existence of regular type N 
solutions with non vanishing twist. 

NOTE ADDED: After this work was completed, we learned of the work by J. 
Bicak and V. Pravda [7], in which the scalar curvature invariants for type N are 
studied. 
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